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Preface 
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and  extremely  helpful  with  his  suggestions  and  questions. 

Most  important , though  I would  like  to  thank  ay  wife,  Letha, 
who  has  encouraged  and  helped  ae  when  she  could  throughout  the 
whole  project,  and  especially  in  the  last  few  hectic  weeks. 


John  T.  Small,  Jr. 


Contents 

Page 

Preface ^ 

List  of  Figures v 

List  of  Tables viH 

Notation ix 

Abstract.  *1 

I.  Introduction 1 

Background  1 

The  Problem 

Purpose  and  Scope  of  the  Study  .....  2 

II.  Theory 4 

Energy  Balance  ^ 

Insolation 5 

Atmospheric  Thermal  Radiation  .....  6 

Thermal  Radiation  of  Objects 8 

Evaporative  Heat  Transfer  ...  9 

Convection . 

Conduction 10 

III.  Models 11 

Sun  Model  11 

Atmospheric  Thermal  Radiation  Model  .........  13 

Air  Temperature  Model  13 

Ground  Model 13 

Tank  Model 17 

Leaf  Model 19 

IV.  Computer  Considerations 23 

Input  Data 23 

Standard  Conditions  24 

Plot  Identification  Numbers  • 25 

V.  Results 26 

Diurnal  Temperature  Curves 26 

Delta-T  Curves.  .......  27 

Parameter  Sensitivity  Curves  and  Analysis  . 28 

VI.  Conclusions  and  Recommendations 55 

iii 


— 


Contents 


2*0. 

Bibliography  37 

Appendix  At  Comparison  of  Ground  Conduction  Model  to 

Analytical  Solution  Using  Sinusoidal  Input.  ....  39 

Appendix  B i Computer  Progress  TEMPS  and  DELTA  .........  63 

Appendix  Ci  Plots  fro*  Prograa  TEMPS  for  all  Parameters  ....  87 

Appendix  Di  Plots  from  Program  DELTA  for  all  Parameters  ....  209 

Appendix  Ei  Sample  Target/Background  Physical  Characteristics  . Z$U 

Vita 257 


( List  of  Figures 

Figure  Page 

1 The  Solar  Spectrum j outside  the  earth's  atmosphere, 

at  14,000  ft.  altitude,  and  at  sea  level 6 

2 Spectral  radiance  of  clear  sky  from 

Elk  Park,  Colorado,  several  elevations  7 

3 Spectral  radiance  of  clear  sky  from 

Cocoa  Beach,  Florida,  several  elevations  7 

4 Variation  of  the  spectral  radiance  of 

the  zenith  sky  with  ambient  air  temperature 8 

5 The  spectral  radiance  of  the  sky  covered  with 

cirrus  clouds  at  several  elevations  8 

6 Spectral  radiance  of  black  body  sources  at 

various  temperatures  as  a function  of  wavelength  ....  9 

7 Energy  balance  terms  for  the  ground  model  14 

8 Verification  of  Eq.  (ll) 1 5 

( 9 The  concept  of  energy  impulse  solution  to  the 

simplified  basic  heat  conduction  equation 17 

10  Energy  balance  terms  for  the  tank  model 18 

11  Energy  balance  terms  for  the  leaf  model 19 

12  Verification  of  Eq.  (17) 21 

13  Explanation  of  Plot  Set  Numbers 25 

14  Sample  Diurnal  Temperature  Curve  for  Summer 

Conditions  at  50°  Latitude  26 

15  Sample  Delta-T  Plot 28 

16  Changes  in  Daytime  Delta-T  caused  by  variations 

in  mean  temperature 30 

17  Changes  in  Nighttime  Delta-T  caused  by  variations 

in  mean  temperature 31 

17(a)  Key  to  Curves  in  Figs.  18(a)  to  1900*  including 

"Standard"  Parameter  Values.  32 

18  Parameter  Analysis)  Daytime 33 

19  Parameter  Analysis)  Nighttime 44 


v 


( 


( 


( 

{ 


Figure 

20  Graph  of  analytical  solution  for  semi-infinite 
slab  heat  conduction  equation  (Hef21»?6).  . . . 

21  Eq.(l4)  solution  to  semi -infinite  slab  problem, 

With  DT-.25  

22  Eq.(l4)  solution  to  semi -inf inite  slab  problem, 

With  DT-.02  

23  TEMPS  solution  with  DT".25 

24  TEMPS  solution  with  DT*.15 


25  Basic  Temperature 

26  Basic  Temperature 

2?  Basic  Temperature 

28  Basic  Temperature 

29  Basic  Temperature 

30  Basic  Temperature 

31  Basic  Temperature 

32  Basic  Temperature 

33  Basic  Temperature 

34  Basic  Temperature 

35  Basic  Temperature 

36  Delta-T  Plots  for 

37  Delta-T  Plots  for 

38  Delta-T  Plots  for 

39  Delta-T  Plots  for 

40  Delta-T  Plots  for 

41  Delta-T  Plots  for 

42  Delta-T  Plots  for 


Plots  of  Mean  Temp 

Plots  of  Absolute  Humidity. 

Plots  of  Wind  

Plots  of  Sun  Strength  . . . 
Plots  of  Ground  Reflectivity 
Plots  of  Ground  Emissivity 
Plots  of  Ground  Diffusivity 
Plots  of  Target  Reflectivity 
Plots  of  Target  Emissivity 
Plots  of  Target  Thickness. 
Plots  of  A.C.  Velocity  . . 

Mean  Temp  

Absolute  Humidity  ..... 

Wind 

Sun  Strength 

Ground  Reflectivity  .... 

Ground  Eaissivity  

Ground  Diffusivity.  .... 


60 


61 

62 

63 

64 
88 
99 

110 

121 

132 

143 

15** 

165 

176 

187 

198 

210 

214 

218 

222 

226 

230 

234 


vi 


I fl I ■-■mp  -■  I.  ..  n.».-  — 

ring*  p&£« 

43  Delta-T  Plots  for  Target  Reflectivity.  .........  238 

44  Delta-T  Plots  for  Target  Eaissivity 242 

45  Delta-T  Plots  for  Target  Thickness  246 

46  Delta-T  Plots  for  A.  C.  Velocity 250 

47  Typical  Reflectance  of  Hater  Surface, 

Snow,  Dry  Soil,  and  Vegetation * . . . . 256 


Table 


List  of  Tables 


ftuce 


I.  Parameters  for  Cloud  Types  12 

II.  Standard  Climatic  Conditions  for  Three 

Model  Locations . . . 24 

III.  Parameter  Identification  Numbers 25 

IT.  Summary  of  Parameter  Strengths  for  Eleven 

Hypothetical,  Geographic,  and  Climatic  Conditions.  . . 55 

V.  Overall  Day  and  Night  Parameter  Strength  Rating.  ...  56 

VI.  Sample  Target/Background  Characteristics  255 


viii 


( 


( 


A (Temp) 

Cd 

C 

P 

CTEMPA 

CTEMP1 

Cv 

^v,in 

C 

v,out 

D 

4T 

d 

z 

E 

Ev 

Enet 

Esurf 

Esat 

es  or  ea 

* 

hc 

^c, force 
^c,free 


Notation 

Portion  of  Grashof  Number  that  represents  buoyancy 
of  fluid 

Heat  Transfer  due  to  conduction 

Thermal  Heat  Capacity  at  constant  pressure  (Cal/gm-°C) 

Computer  variable  for  air  temperature 

Component  used  in  calculating  the  air  temperature 

Heat  Transfer  due  to  convection 

Heat  convection  inside  the  tank  model 

Heat  Convection  outside  the  tank  model 

Diffusivity  (cm  /min) 

Temperature  difference  between  two  objects 

Ratio  of  molecular  weights  also  sun  declination  angle 

Enissivity,  or  small  time  increment 

Net  energy  deposited  on  a surface  (cal/cm^-min) 

Total  energy  in  a volume  V of  substance 
Net  energy  stored  in  an  impulse  in  time  dt 
Energy  per  unit  volume  at  the  surface  of  the  ground 
Saturated  water  vapor  pressure  (millibars) 

Vapor  pressure  at  earth's  surface  or  in  the  air 

(respectively) 

Latitude  in  degrees 

Viscocity  of  the  air 

Heat  transfer  coefficient  (general) 

Heat  transfer  coefficient  due  to  forced  convection 
Heat  transfer  coefficient  due  to  free  convection 


ix 


K 

°K 

LE 

Pr 

R 


RH 

Re 


RSOIL 

S 


S 

P 


P 

P, 

e 

a 

u 

X 

Za 


VonKarmon’s  constant  (.41 ) 

Degrees  Kelvin 

Energy  loss  due  to  evaporation  of  moisture 
Prandtl  number  (Approximately  ■ .72  for  air) 

Total  radiation  of  a body 

Thermal  radiation  gain  due  to  surrounding  radiating 
objects 

Thermal  radiation  loss  of  objects 

Relative  Humidity  (0  to  100%) 

Reynolds  number 

Reflectivity  of  the  soil 

Insolation  due  to  Sun  at  earth’s  surface 
corrected  for  clouds  and  orientation  angle 

Atmospheric  Thermal  radiation 

Total  sun's  energy  on  an  object  which  is  perpen- 
dicular to  the  sun's  rays 

Insolation  on  an  object  oriented  parallel  to  the 
ground  (uncorrected  for  clouds) 

Density  of  an  object 

Density  of  the  air 

Time 

Stefan-Boltzmann  constant 
velocity  of  the  air  (cm/sec) 
distance  in  cm 

Sensor  height  (for  measuring  wind  or  tempera ture ) 
Roughness  parameter  (about  .025  cm  for  smooth  surfaces) 


x 


GXP/PH/77-12 


An  investigation  of  the  thermal  signatures  of  natural  objects  la 
undertaken.  Using  the  principle  of  energy  balance,  diurnal  temperature 
models  are  developed  for  the  ground,  a vehicle  (such  as  a tank),  and 
the  leaves  of  a tree.  These  nodela  are  included  in  a computer  program 
designed  to  simulate  daily  variations  in  conductive,  radiative,  and 
convective  heat  transfer  processes. 

Environmental  conditions  are  changed  by  altering  the  program  ln- 

j 

puts,  which  include | latitude,  day,  wind  speed,  mean  air  temperature, 
cloud  type,  total  insolation,  absolute  and  relative  humidity,  atmospheric 
pressure,  and  particulate  concentration.  The  program  is  iterated  for 

several  values  of  a single  parameter  such  as  one  of  the  environmental 

« 

inputs,  or  one  of  21  physical  characteristics  of  the  models  (for  example} 
ground  reflectivity,  tank  thickness,  or  leaf  transpiration  rate).  There 
are  two  types  of  output  from  the  program}  simultaneous  plots  of  the  temp- 
eratures of  all  the  models  for  a 30-hour  period,  and  curves  representing 
24-hour  periods  of  the  temperature  difference  (Delta-T)  between  the  tank 
and  the  ground.  The  effects  of  changing  a parameter  are  analysed  by 
comparing  the  various  Delta-T  curves. 

In  this  study,  11  scenarios  (representing  three  geographic  locations 
and  a variety  of  climatic  conditions)  are  used  to  examine  the  following 
parameters}  mean  temperature,  absolute  humidity,  wind  speed,  cloud  oover, 
ground/ target  reflectivity  and  enlsslvity,  ground  diffuslvity,  target  thick- 
ness, and  tank  air  conditioning  velocity.  Analysis  of  the  Delta-T  curves 
indicates  that  overall,  the  ground  enisaivlty  appears  to  affect  the  temp- 
eratures the  most,  while  the  amount  of  insolation  affects  them  the  least. 
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A THEORETICAL  ANALYSIS  OF  CHANGES 


I 


IN  THERMAL  SIGNATURES  CAUSED  BY 
PHYSICAL  AND  CLIMATOLOGICAL  FACTORS 

I.  Introduction 


Background 

Although  thermal  imaging  systems  have  been  used  since  the  1940’s, 
accurate  nays  of  predicting  the  thermal  signatures  of  target  backgrounds 
have  not  been  developed  (Ref  1 1267).  Modem  Infrared  imaging  devices 
oan  detect  the  temperature  differences  between  objects  such  as  trees, 
fields,  roads , tanks  and  trucks,  and  can  convert  these  thermal  contrast 
differences  into  a visual  display  on  a screen.  Because  such  equipment 
is  an  integral  part  in  the  development  of  new  missle  systems  like  the 
Air  Force  Imaging  Infrared  Maverick  and  the  Amy  Hellflre,  it  is  impor- 
tant to  have  an  understanding  of  what  environmental  factors  drive  the 
temperature  differences  of  the  objects  detected  by  these  infrared  de- 
vices. A knowledge  of  the  physics  behind  thermal  signatures  can  increase 
the  insight  into  what  infrared  images  should  look  like  under  various  con- 
ditions. This  Insight  should  be  helpful  not  only  to  scientists  gathering 
data  on  infrared  image  characteristics,  but  also  to  system  operators 
trying  to  identify  particular  targets. 


i 
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The  Problem 

In  spite  of  the  accuracy  of  current  infrared  systems,  incorrect 
target  signature  prediction  has  still  caused  problems.  For  example, 
in  one  test  of  the  Maverick  missle,  inaccurate  aiming  occured,  even 
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though  the  crew-members  had  studied  videotapes  of  Infrared  Images  made 
In  Western  Europe.  These  tests  were  sonltored  by  the  Office  of  the  Assis- 
tant Chief  of  Staff  for  Studies  and  Analysis,  Air  Force  Headquarters, 
Washington,  D.C. , After  reviewing  the  results,  a Member  of  that  office, 

Lt  Col  Carl  Case,  suggested  that  a study  be  started  to  Investigate  the 
physical  processes  of  thermal  heating,  with  the  premise  that  a better 
understanding  of  the  physics  Involved  may  lead  to  better  prediction  of 
target  signature  characteristics  under  varying  environmental  conditions. 

Unfortunately , current  work  on  thermal  signatures  has  been  re- 
stricted to  very  specific  areas,  and  has  been  concentrated  on  a single 
Item  such  as  a tank  (Ref  2),  or  the  ground  (Ref  3«  Ref  4).  There  does 
not  seen  to  be  a study  using  several  items  In  a scenario  to  determine 
how  their  interrelated  temperature  differences  would  change  under  a 
variety  of  environments. 

Purpose  and  Scope  of  the  Study 

This  analysis  applies  heat  transfer  principles  in  a computer  model 
to  determine  the  key  driving  factors  behind  the  daily  and  seasonal  temp- 
erature differences  of  various  objects  and  their  backgrounds.  The  models 
that  are  developed  simulate  heat  transfer  processes  at  the  surfaces  of 
the  ground,  a tank,  and  the  leaf  of  a deciduous  tree.  The  study  Is  kept 
very  general,  since  specific  results  have  only  limited  applicability 
and  nay  require  lengthy,  detailed,  computer  coding.  By  relying  on  a 
theoretical,  physical  appraisal,  as  opposed  to  a complex  mathematical  sol- 
ution, the  output  of  the  analysis  will  be  a series  of  graphs  showing  how  ther- 
mal signatures  vary  with  certain  general  parameters  (to  be  described  later) 
rather  than  a precise  method  of  predicting  temperatures  at  speoifio  times. 
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Two  computer  programs  will  be  used.  The  first  program  (TEMPS) 
calculates  and  graphs  a 30  hour,  diurnal  temperature  profile  of  the 
air  temperature,  and  the  tank,  ground,  and  tree  models.  The  second 
program  (DELTA)  computes  the  same  temperature  cycles  and  plots  out 
the  difference  between  the  tank  and  the  ground  temperature  from  the 
6th  to  the  30th  hour. 

Both  programs  are  set  up  to  iterate  four  times  in  each  run, 
and  a selected  parameter  is  varied  for  each  iteration.  The  tank/ground 
temperature  difference  (AT  ) is  noted  at  the  14th  and  24th  hour 
(representing  times  of  high  and  low  heat-loading,  respectively)  for 
all  values  of  the  various  parameters.  These  data  are  then  graphically 
compared  and  analyzed  to  determine  the  relative  effect  that  each 
parameter  has  on  A T . 
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II.  Theory 

The  tenperature  of  an  object  arises  primarily  due  to  the  kinetic 
energy  of  its  molecules.  In  a natural  setting,  this  energy  is  dis- 
tributed through  a complex  energy  balance  with  the  environment.  This 
balance  changes  continually,  depending  upon  the  object  Itself,  its  lo- 
cation, the  time  of  year,  the  time  of  day,  and  the  existing  meteorolo- 
gical conditions. 

Energy  Balance 

If  the  amount  of  random , or  thermal  energy  per  volume  of  sub- 
stance can  be  determined,  then  the  temperature  of  the  object  can  be 
calculated  through  the  relation i 

T = h 

PCpV  (1) 

where  Ev  is  the  thermal  energy  contained  in  a volume  V of  material 
with  density  p and  a thermal  heat  capacity  at  constant  pressure,  Cp  . 
Because  of  the  complex  interrelationships  of  internal  and  external 
factors  it  is  extremely  difficult  to  determine  the  precise  amount 
of  energy  flow  into  and  out  of  an  object  at  a specific  time.  As  a 
result,  determining  exact  temperatures  from  the  net  energy  remaining 
with  the  object  is  also  difficult. 

However,  by  setting  up  an  equation  which  represents  the  energy 
flow  of  various  heat  transfer  processes  (Ref  50),  a model  can  be 
generated  to  determine  the  relative  significance  of  those  processes. 

The  above  technique  is  used  in  the  models  of  this  study  and  the  fol- 

k 

“ _ ' 
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lowing  energy  balance  equation  is  assumed i 

E = $ + VW  LE±  Cv±  Cd  (2) 

a 

(All  terms  in  gram-calories  per  cnr  per  min) 

Where  E « the  net  energy  being  deposited  at  the  surface 

S » the  solar  insolation  (direct,  scattered  and  reflected) 

* incident  thermal  radiation  from  the  atmosphere 
Rj  - thermal  radiation  gain  due  to  surrounding  objects 
R0«  thermal  radiation  loss  of  the  object 

IE  ■ energy  loss  due  to  evaporation  of  moisture  from  surface 
Cv  ■ energy  loss  (gain)  due  to  free  or  forced  convection 
Cj  " energy  loss  (gain)  due  to  thermal  conduction. 

In  all  cases  the  convention  of  a plus  sign  means  flow  into  the  body. 

Other  terms  that  represent  such  processes  as  heat  generated  by 
internal  sources  (a  running  engine,  for  example)  can  be  added  to  Eq  (2), 
but  will  not  be  lnoluded  explicitly  in  the  models  of  this  study.  (An 
"engine"  can  be  simulated  by  assuming  a higher  interior  temperature). 
Because  of  various  conditions,  some  terms  in  Eq  (2)  may  not  be  used 
in  certain  models.  Descriptions  of  the  energy  balance  terms  and  their 
associated  heat  transfer  processes  follows. 

insolation  (S) 

The  main  thermal  energy  supply  to  objects  on  the  earth  comes 
from  the  sun.  The  sun’s  radiation  can  be  approximated  by  a 6000°  K 
black  body,  and  at  the  outer  edge  of  the  earth's  atmosphere  it  has  an 
intensity  of  about  2.0  cal/cm2-min  (Ref  6).  After  propagating  through 
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Fig.  1.  The  solar  spectrum j outside  the  earth’s  atmosphere, 
at  14,000  ft  altitude,  and  at  sea  level  (Ref  5*^9) 

the  atmosphere,  however,  absorption  and  scattering  has  reduced  it 
by  one  half,  and  the  spectral  characteristics  are  modified  consider- 
ably. Figure  1 shows  that  the  largest  part  of  incident  energy  is  in 
the  visible  and  near  infrared  region,  thus,  for  radiation  heat  pur- 
poses, the  spectral  properties  of  objects  at  these  wavelengths  are 
critical. 

The  amount  of  scattering  and  absorption  is  further  related  to 
the  distance  of  transmission  through  the  atmosphere  and  this  depends 
on  the  time  of  year  and  time  of  day  for  each  point  on  the  globe.  Hence, 
to  arrive  at  the  amount  of  solar  energy  impinging  upon  an  object  at 
a certain  time,  some  of  the  factors  that  must  be  known  arej  latitude, 
astronomical  (or  sun)  time,  amount  of  water  vapor  in  the  atmosphere 
( a key  absorber  in  the  infrared),  amount  of  particulate  matter  in  the 
air  (a  cause  of  scattering)  and  the  type  and  amount  of  cloud  cover. 

Atmospheric  Thermal  Radiation  ( St ) 

The  atmosphere  absorbs  and  scatters  the  sunlight  which  passes 
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Fig.  2.  Spectral  radiance  of 
clear  sky  from  Elk  Park,  Col- 
orado, several  elevations. 
(Ref  7*1315) 


MCRONS  WAVELENGTH 

Fig.  3.  Spectral  radiance  of 
clear  sky  from  Cocoa  Beach,  Flo 
rida,  several  elevations. 

(Ref  ?i1316) 


through  it.  Additionally  the  molecules  re-radiate  thermally  at  long 
wavelengths  (8-12  microns)  like  a grey  body  with  an  eaissivity  that 
depends  on  the  amount  of  gasses  (particularly  water  vapor  and  carbon 
dioxide )present.  The  sky  radiance  is  highly  dependent  upon  the  el- 
evation angle,  since  the  effective  thickness  of  the  atmosphere  is 
least  when  measured  directly  overhead.  The  influence  of  the  climate 
and  elevation  angle  is  shown  in  Figs.  2 and  3,  where  the  drier  air 
in  Fig.  2 exhibits  a lower  peak.  Since  the  molecular  activity  depends 
greatly  upon  the  temperature,  the  ambient  air  temperature  also  has  a 
strong  influence  on  the  radiance,  as  shown  in  Fig.  4.  The  presence 
of  clouds  (Fig.  5)  primarily  tends  to  increase  the  amount  of  radiance 
in  the  zenith  direction,  which  also  increases  the  total  amount  of  sky 
thermal  radiation  received  by  the  earth. 
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Fig.  4.  Variation  of  the  spectral  Fig.  5.  The  spectral  radiance  of 
radiance  of  the  zenith  sky  with  the  sky  covered  with  cirrus  clouds 

ambient  air  temperature  (fief  7 *1316)  at  several  elevations,  (fief  7*1318) 


Thermal  Radiation  of  Objects  ( Rj  and  ftp  ) 

All  objects  above  absolute  zero  give  off  thermal  radiation  due 
to  Random,  molecular  activity.  The  total  amount  of  radiation  given 
off  is  determined  by  the  Stefan-Boltzmann  law 

* - (3) 

Where  e is  the  emisslvity  of  the  object  which  characterizes  its 
emission  compared  to  that  of  a perfect  black  body,  <r  is  the  Stefan- 
Boltsmann  constant  of  7.93  x 10”^  cal  per  cm2  per  °K^per  minute 
(5.67  x 10"  org  per  CIS  t>er  °K  per  second)  and  T is  the  temperature 
of  the  object  in  °K.  The  spectral  distribution  of  thermal  radiation  is 
governed  by  the  Planck  radiation  law  and  the  wavelenth  of  peak  emmisslon 
for  most  natural  objects  is  in  the  8 to  10  micron  region  as  shown  in 
Fig.  6.  The  constant  loss  of  energy  by  thermal  radiation  is  one  of 
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WAVELENGTH  IN  MICRONS 

Fig.  6.  Spectral  radiance  of  black  body  sources 
at  various  temperatures  as  a function  of  wave- 
length (Ref  3*62). 

the  major,  natural  mechanisms  for  temperature  regulation  in  objects. 

R0  represents  the  outgoing  thermal  radiation  of  the  object  con- 
sidered and  R j is  the  incoming  thermal  radiation  due  to  other  near- 
by thermal  sources  other  than  the  atmosphere. 

Evaporative  Heat  Transfer  ( It  ) 

When  one  gram  of  witer  is  converted  from  liquid  to  vapor,  approx- 
imately 600  calories  of  heat  is  needed.  Thus,  the  evaporation  pro- 
cess that  occurs  at  the  surface  of  a leaf,  or  at  the  top  of  the  soil 
is  capable  of  removing  a significant  amount  of  thermal  energy,  and 
reducing  the  temperature  at  the  surface.  The  amount  of  evaporation  is 
governed  by  natural  factors  such  as  the  diffusion  resistance.  This  term 
characterises  the  rate  of  fluid  transport  to  the  surfaoe,  and  since  its 
value  can  vary  greatly,  so  can  the  quantity  of  evaporative  cooling. 

Convection  ( ) 

Convection  is  a method  of  heat  transfer  caused  by  the  mixing  of 
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fluids  of  different  temperatures  (and  hence  different  densities). 
During  convection,  the  flow  of  heat  across  a temperature  gradient 
is  generally  enhanced  because  the  fluids  on  each  6ide  of  the  grad- 
ient are  nixed  in  the  process.  There  are  two  types  of  convection, 
free  convection,  where  the  movement  is  caused  by  the  temperature -re- 
lated buoyancy  of  the  air,  and  forced  convection,  where  a large-scale 
mass  movement,  such  as  wind,  is  involved.  Both  types  of  convective 
heat  transfer  depend  on  the  shape  and  orientation  of  the  object,  as 
well  as  the  temperature  difference  between  it  and  the  surroundings. 


Conduction  (Cd) 

Heat  transfer  by  conduction  depends  on  the  direct  exchange  of 
kinetic  energy  between  molecules.  The  rate  of  this  transfer  is 
highly  dependent  upon  the  physical  makeup  (thermal  conductivity) 
of  the  substance  and  the  temperature  difference  across  the  conduct- 
ing path.  Hence  the  contribution  of  conduction  to  the  overall  energy 
balance  may  vary  widely. 

It  is  assumed  that  the  time  dependence  of  conductive  heat  trans- 
fer follows  the  classic  heat  equation  (Ref  81III-I8) 


i T 
d 6 


= D 


a'T 

ax* 


(4) 


Where  T is  the  temperature,  X is  the  distance  below  the  surface, 

0 is  the  time,  and  D is  the  thermal  diffusivity. 

Simplifying  the  geometry  of  the  models  will  allow  the  use  of  a 
one  dimensional  solution  to  Eq  (4).  Thin,  highly  conductive  objects 
will  reach  equilibrium  quickly  so  Sq  (4)  will  not  be  necessary.  In  such 
cases  temperature  can  be  determined  by  energy  balance  alone. 
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III.  Models 

In  order  to  analyze  the  effect  that  various  physical  parameters 
and  climatological  conditions  have  on  Delta-T  (the  temperature  differ- 
ence of  am  object  relative  to  its  background)  it  is  necessary  to  develop 
models  for  at  least  two  different  media.  Because  of  the  military  inter- 
est in  using  infrared  systems  for  tank  detection,  the  models  used  are  in- 
tended to  give  a rough  representation  of  tank  thermal  signatures.  These 
models  encompass)  the  ground,  a tank,  and  the  leaves  of  a deciduous 
tree.  Each  of  these  will  be  discussed  in  this  section)  however,  first  will 
come  a description  of  some  other  models  needed  to  simulate  the  sun's 
strength,  the  atmospheric  thermal  radiation,  and  the  air  temperature. 

Sun  Model 

The  solar  zenith  angle,  with  respect  to  the  normal  at  the  earth's 
surface,  i , can  be  calculated  from  the  following  equation  (Ref  3»58)» 

Cos(z)  = Sin(0)Sin(d)+Cos(0)Cos(d)Co$(h)  (5) 

Where  0 is  the  latitude  ( + is  north),  d is  the  sun's  declination, 
and  h is  the  sun's  hour  angle.  The  angle  0 is  given  as  one  of 
the  inputs  to  the  computer  program,  and  9 is  calculated  from 
a curve  fit  of  the  data  in  the  Air  Almanac  (Ref  9)  using t 

d - [jg  x 52.88J  - 26.44  (fort  0<Day<l83) 
dm  26.44  - jjjg  x 52.88]  (fort  184  < Day  <365) 

Where  Day  is  the  Julian  date. 


(6a) 

(6b) 
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The  hour  angle  is  assumed  to  be  set  so  that  12i00  noon  occurs  when 
h “ 90°.  Thus,  there  is  no  correction  for  local  time  (longitude). 


The  strength  of  the  sun  on  a flat  surface  at  the  ground  and 
perpendicular  to  the  sun's  rays  is  calculated  using  an  equation  de- 
veloped after  Brooks  (Ref  10*264) t 


2 

cal  /cm  -min 


Where  p “ air  pressure  (mb) 

w ■ total  precipitable  water  vapor  in  the  atmosphere  in  the 
zenith  direction  (mm) 

d «*  concentration  of  particulate  matter  (particles/ca^) 
m « secant  of  the  zenith  angle. 

To  correct  this  value  for  a horizontal  surface,  Sp  is  multi- 
plied by  the  result  of  Sq.  (5);  $0  * Sp  Cos  (z),  A final  adjustment 

for  the  amount  of  insolation  comes  from  the  effect  of  clouds,  and 
is  contained  in  an  empirical  formula  by  Haurwitz  (Ref  11 tll2 ) t 

S " 8 0^4]  6XP  [ "® (b~* 059 }]  (8 } 

Where  m is  the  secant  of  the  zenith  angle,  a,  and  b are  parameters 
based  on  cloud  type  (See  Table  i). 


Table  I 

Parameters  for  Cloud  Types  (Ci,  cirrus j Cs,  cirrostratusj  Ac,  altocu- 
mulus* As,  altostratusi  Sc,  stratocumulusi  St,  stratus*  Ns.  nlmbostratus  ) 


Ci 

Cs 

Ac 

A8 

Sc 

St 

Ns 

Fog 

a 

82.2 

87.1 

52.5 

39.0 

34.7 

23.8 

11.2 

15.4 

b 

.079 

.148 

.112 

.063 

.104 

. .Iff. 

-.167 

.028 

(From  Ref  3 >56) 
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Atmospheric  Thermal  Radiation  Model 


As  previously  mentioned,  there  are  many  variables  affecting  the 
thermal  radiation  of  the  atmosphere.  Complex,  analytical  solutions 
have  been  developed  for  this  problem,  but  they  require  detailed  data 
for  the  calculations.  For  this  study,  a simple  empirical  formula, 
developed  by  Angstrom  for  a clear  sky  is  used  (Ref  12 )i 

S,  - o-T4  (0.81-.2411  x 10_,°52  **)  (9) 

Where  ea  is  the  partial  pressure  of  water  vapor  at  ground  level  in 
millibars  (calculation  of  ea  will  be  discussed  later). 

Air  Temperature  Model 

An  exact  model  for  dally  air  temperatures  is  difficult  due  to 
the  effects  of  such  meteorological  conditions  as  winds  and  fronts. 
Nevertheless  an  arbitrary  diurnal  cycle  was  chosen  using  the  temp- 
erature profiles  noticed  in  the  literature  (Ref  13i?8j  Ref  14j109» 

Ref  15i43-46).  Care  was  taken  to  attempt  to  preserve  the  times  of 
occurrence  of  maximum  and  minimum  temperatures. 

CTEMPA  - CTEMP1  + 10  Sin(  n (1 .05+  ^ ))  + 2.5  Sin(  n (l  .5+  ^))  (10) 

Where  CTEMPA  is  the  air  temperature  in  °C  at  time,  t,  CTEMP1  is  a 
"DC"  term  to  shift  the  temperature  according  to  the  season  and  lat- 
itude, and  t is  the  time  of  day, in  hours, since  midnight. 

Ground  Model 

The  ground  is  assumed  to  be  a barren,  semi-infinite  slab,  composed 
of  thermally  homogeneous  and  isotropic  soil.  This  simplifies  the 
ground  heat  conduction  solution,  so  that,  providing  the  energy  balance 
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Fig.  7.  Energy  balance 
terms  for  the  ground  model. 

remains  uniform  along  the  boundary,  the  direction  of  heat  flow  will 
be  one  dimensional.  Figure  7 identifies  the  energy  balance  terms 
that  are  assumed  to  be  significant  for  the  ground  model. 

The  insolation  ( S ) and  atmospheric  long-wave  radiation  ( Sj  ) 
are  equal  to  the  values  computed  from  their  respective  models,  depend- 
ing on  the  time  of  day  and  conditions.  The  rate  of  thermal  heat 
loss  ( R0  ) from  the  surface  is  based  upon  Eq.  (3),  using  the  ground 
temperature  calculated  in  the  previous  program  time  Increment  and  an 
emissivity  that  is  set  at  the  beginning  of  the  computer  run. 

Use  of  analytical  convective  heat  transfer  equations  for  the 
ground  model  fail  due  to  a dependence  upon  the  "characteristic  length" 
(the  dimension  of  the  object  in  the  direction  of  flowi  which  is  infinite 
for  the  ground  model).  Van  Bavel  (Ref  17(85)  was  also  unable  to  derive 
an  applicable  ground  model  formula.  Thus,  the  forced  conveotion  coef- 

, fieients  in  this  study  axe  from  an  empirical  formula,  Eq  (11),  based 

t 

on  convection  data  over  a flat  surface  (Ref  18i538).  The  free  convec- 
tion coefficient  also  comes  from  Eq  (ll ),  with  the  wind  speed  set  at 
sero.  This  value  was  checked  by  Inserting  it  in  a standard  equation 
(Ref  16»296)  and  solving  for  the  characteristic  length.  The  resulting 
dimension  of  85  cm  is  large  compared  to  typloal  data  (Ref  16)  and  so  for 
the  "infinite"  ground  it  appears  that  the  empirical  formula  is  also 
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a reasonable  approximation  for  free  convection.  The  formula  used  isi 

S - «•*>  - f2^]  ’<•■»]  >*  ■*--’>  O'  > 

Cal/cm2-min-°C 

Where  Speed  is  the  velocity  of  the  air  over  the  surface. 

The  expression  for  (if  ) is  used  from  Van  Bavel  (Ref  I?i7) 

2 

IE  - (L)  -fi-1  k,“  (6o(es  - ea))  (12) 

P [ln(za/Zo )f 

Cal/cm  -min 

Where  L “ Latent  heat  of  Vaporization  (*  600  Gal/gm) 

■ Density  of  air  (gra/cm-^) 
k ■ Von  Karman's  Constant  (.41) 

& - Water/air  molecular  weight  ratio  (.622) 
u « Velocity  of  air  (cm/sec) 

P ■ Atmospheric  pressure  (mb) 

Za“  Air  sensor  height  (100  cm) 

ZQ-  Roughness  length  (.025  cm) 
es/ea  ■=  Water  vapor  pressure  at  gound/alr  (mb) 

Za  was  arbitrarily  chosen  to  be  100  cm  as  a representative  height  at 
which  temperature  can  be  measured,  the  value  for  Z0  is  from  Van  Bavel 
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(Ref  I?i85)  and.  Is  representative  of  bare  soil.  Saturated  vapor 
pressure  (esat)  is  calculated  from  an  empirical  formula  by  Brunt 
(Ref  19«104)i 


esat  " 


(5423) 


(eo) 


(mb) 


(13) 


Where  Tg^  is  the  temperature  of  the  saturated  vapor  in  °K  and  eo  is 
the  water  vapor  pressure  of  saturated  air  at  0°C  ( 6.1071  mb).  This 
vapor  pressure  must  then  be  corrected  for  relative  humidity  using  the 
following  (Ref  20il62)i 


ea  ■ esat(RH/lOO) 


(1*0 


Where  RH  is  the  relative  humidity  from  0 to  100%.  An  expression  sim- 
ilar to  Sq  (14)  is  used  for  es  except  that  instead  of  relative  hum- 
idity the  moisture  content  of  the  ground  (in  percent  by  weight)  is 
substituted  for  RH. 

At  this  point  in  the  model  all  the  energy  terms  for  the  ground 
are  summed  ( S + S|  - R -IE  ± C¥)  to  give  a net  flow  of  energy  (Enet  in 
cal/cm^-min)  which  is  then  used  to  calculate  the  amount  of  thermal 
heat  conduction.  Carslaw  and  Jaeger  (Ref  21 ) have  developed  many 
analytical  solutions  to  heat  conduction | however,  most  of  the  input 
heat  fluxes  are  explicit  functions  of  time.  Because  of  the  activity 
of  clouds  and  winds,  the  net  heat  flux  into  the  ground  may  not  be  a 
smooth,  continuous,  curve,  so  an  approximation  to  the  solution  can 
be  assumed  as  follows  1 

Bq  (4)  becomes  1 _L_  « -1—  D 

6 X 

Yielding  1 X (13) 
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Pig.  9.  The  concept  of  the  energy  impulse  solution  to  the  sim- 
plified basic  heat  conduction  equation. 

Eq.  (13)  represents  an  approximate  distance,  X that  a heat  wave  will 

travel  through  a medium  with  diffusivity,  0 in  a time  6 . 

If  the  net  energy  flow,  Enet  is  "stored  up"  and  deposited  as 

an  impulse  every  "dt"  minutes,  then  each  impulse  will  diffuse  into 

the  soil  between  intervals  (see  Fig.  9).  This  will  cause  a net 

"piling -up"  of  energy  at  the  surface.  The  net  amount  of  energy  that 

remains  at  the  surface  after  N impulses  can  be  expressed  as: 

Esurf  ■ E (*<■?]  ♦pMt(dt)  ] (14) 

m“l  y 0 (m  (dt)) 

This  net  energy  gain  is  then  used  in  Eq  (l ) to  calculate  the  net 
temperature  difference  at  the  earth’s  surface.  Eq.  (14)  has  been 
tested  against  the  analytical  solution  for  a sinusoidal  heat  flux 
input  (Ref  21«?6)  with  good  results  (see  Appendix  A).  The  computer 
ground  model  has  one  additional  computation  included  to  allow  for 
the  reverse  flow  of  energy  when  the  sub-terrainian  temperature  is 
greater  than  at  the  surface. 

Tank  Model 

The  tank  model  consists  of  a hollow,  rectangular  box  assumed 

to  be  made  of  Iron,  3 in  (7.^2  cm)  thick.  The  height  of  the  box  is 
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considerably  smaller  than  the  length  or  width  so  that  the  heat  flux 
contribution  from  the  sides  is  considered  to  be  negligible.  The  top 
of  the  box  facing  the  sky  is  the  surface  of  primary  interest  for 
temperature  considerations.  Although  the  bottom  surface  receives 
long-wave  radiation  from  the  ground,  it  is  assumed  that  this  heats  the 
lower  surface,  causing  it  to  radiate,  but  only  enough  to  approximately 
balance  the  radiation  from  the  heated  upper  surface.  (Under  high 
heat -loading  this  simplification  may  cause  the  computed  temperature 
of  the  upper  surface  to  be  erroneously  high).  The  dominant  heat  factors 
considered  are  shown  in  Fig.  10. 

The  thickness  of  the  skin  is  such  that  it  reaches  equilibrium 
within  the  time  increments  of  the  overall  program  (usually  15  minutes). 
This  can  be  verified  by  solving  Eq  (13)  for  0 and  using  the  diffusiv- 
ity  of  iron  and  substituting  X=  7.62  cm.  In  this  case,  0 turns  out 
to  be  about  8 minutes  for  the  heat  wave  to  travel  through  the  skin. 

The  computations  for  S , Sj  , R0  and  CY  0U|  are  the  same  as  for 
the  ground  model.  The  addition  of  the  term  Cyjn  is  to  include  the 
simulation  of  heat  exchange  with  the  atmosphere  inside  the  tank.  The 


Fig.  10.  Energy  balance 
terms  for  the  tank  model. 
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atmosphere  is  assumed,  to  be  air  conditioned  at  2?°C,  and  is  kept 
moving  at  a velocity  of  44?  cm/sec.  Although  this  model  appears  to 
be  extremely  simplified,  It  roughly  parallels  a very  successful  tech- 
nique used  by  the  Any  Night  Vision  Laboratories,  where  the  temp- 
erature signatures  of  tanks  and  trucks  are  area-weighted  to  arrive 
at  a single , overall  A T that  is  then  used  to  characterize 
each  vehicle  (Ref  22).  Computation  of  Cvjn  is  the  same  as  forCv0ut  , 
using  Eq.  (ll). 

Leaf  Model 

The  leaf  model  is  assumed  to  be  a flat  plate  that  is  suspended 
parallel  to  the  ground  and  is  fully  exposed  to  the  sky  and  ground  for 
the  complete  30  hour  cycle.  This  would  be  typical  of  an  outer  leaf 
located  on  a lower  branch  of  a tree.  The  thickness  is  small,  so  that 
a time  independent  solution  is  valid.  Figure  11  indicates  what  are 
considered  to  be  the  primary  heat  transfer  terms  for  the  leaf.  Note 
that  both  sides  make  significant  contributions  to  the  energy  balance 
equation. 


Fig.  11.  Energy  balance  terms 
for  the  leaf  model. 
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The  upper  surface  of  the  leaf  receives  the  direct  portion  of 
insolation  ( $ ) and  the  atmospheric  thermal  radiation  ( Sj  ).  The 
lower  surface  intercepts  not  only  the  thermal  radiation  of  the  ground 
( R { ) i but  also  a portion  of  the  insolation  that  has  been  reflect- 
ed off  the  ground  ( S x Reflectivity  of  the  Ground). 

Outgoing  energy  from  the  leaf's  top  and  bottom  sides  leaves  in 
the  form  of  thermal  radiation  ( Ro),  evaporation  through  transpor- 
ation  (LE)  and  convection  (Cv).  The  processes  of  all  the  heat  ex- 
change terms  are  the  same  as  in  the  previous  models,  except  for  small 
differences  in  the  convection  and  evaporation  terms. 

The  geometry  of  a leaf  is  such  that  it  applies  to  the  basic 
equations  for  heat  transfer  by  convection.  Thus  a distinction  must 
be  made  between  the  free  and  forced  convection  regimes.  The  heat  trans- 
fer coefficient  for  free  convection  over  a flat  plate  is  (Ref  51 100) 


hCffree  - 92^1  (Gr  (Pr)  )*25  (15) 

Where  L is  the  characteristic  length  of  the  plate  (leaf),  Pr  is  the 
Prandtl  number  of  the  air  ( « .72  for  temperatures  between  0 and  100  °C), 
and  Gr  is  the  Grashof  number. 

But  Gr  - A(Temp)  (AT)(L)^  (16) 

Where  A (Temp)  is  a buoyancy  term  that  is  a function  of  the  temperature 
and  pressure.  A curve-fit  relation  for  the  values  of  A(Temp)  has 
been  included  in  the  program  using  data  from  Kreith  (Ref  l6»595)* 


A(Temp)“ 
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Agreement  for  Eq.  (l?)  is  shown  In  Fig.  12.  Using  the  results  of  Eq.  (17) 
leaves  the  following  relationship  for  free  convectionj 


c,free 


Const 


A (Temp)  ( AT  )(Pr) 


.25 


(18) 


For  forced  convection,  the  heat  transfer  coefficient  is  related 
by  (Ref  16  *295 ) s 


c,  force 


Const  (R®)*5(Pr)*33 
L 


Where  Re  is  th3  Reynolds  number  and  is  computed  from; 


(19) 


Re  =•  L (20) 

H- 

Where  p.  is  the  absolute  viscocity  of  the  air  (gm/cm-sec)  and  u is 
the  air  velocity  (cm/sec).  Using  the  ideal  gas  law  to  calculate 
Pj  and  a curve-fit  equation  for  h , fitted  to  within  ,5/S  of  a set  of 
tables  (Ref  18 til)  for  temperatures  between  240  and  330  °K,  the 
forced  convection  coefficient  is  then  computed  £>nce  the  air  temp- 
erature, air  velocity  and  leaf  dimension  are  known)usingj 


h - (Constant)  (Pr),33[~u  Pa  1 *5 

c .force  ^ L 


(21) 
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The  program  compares  the  values  of  hc  free  and  hQ  fQTce  and  the  higher 
of  the  two  is  used  for  computation. 

Evaporation  (IE)  at  the  leaf  surface  is  accomplished  through 

transpiration,  which  depends  upon  the  quantity  of  water  pumped  to 
the  leaf's  surface,  and  is  controlled  by  the  plant's  physiology. 

In  this  process,  moisture  reaches  the  outside  of  the  leaf  through 
pore-like  cells,  called  stomata.  "Guard  cells",  at  the  opening  of 
these  pores , act  as  valves  to  regulate  how  fast  the  water  reaches 
the  surface,  where  evaporation  (and  hence  cooling)  can  take  plane. 

The  amount  of  water  transported  can  vary  from  0 to  0.001  gn/cm  -min 
(Ref  25)994),  depending  upon  the  type  of  plant,  its  "health",  and  the 
heat-loading  conditions.  The  model  in  this  study  is  designed  to 
follow  the  leaf  temperature  cycles  Indicated  by  Gates  (Ref  24i82). 

The  transpiration  rate  is  very  near  aero  at  night,  and  if  winds 
are  light,  the  leaf  can  cool  down  below  the  air  temperature  by  rad- 
iation alone.  After  sunrise,  the  heat-loading  increases,  as  does 
the  leaf  temperature,  and  transpiration  starts  when  the  leaf  temp- 
erature reaches  the  air  temperature.  This  condition  is  maintained 
by  regulating  the  transpiration  rate  to  provide  sufficient  cooling 
moisture.  If  the  heat-load  demand  exceeds  a preset  maximum  trans- 
piration rate,  the  leaf  temperature  will  rise  above  the  air  temper- 
ature to  an  equilibrium  based  on  energy  balance.  Later  in  the  day, 
transpiration  decreases  as  the  heat-loading  goes  down  and  the  leaf 
temperature  reaches  the  air  temperature.  Eventually,  at  night,  the 
transpiration  rate  drops  to  zero,  and  the  leaf  temperature  is  again 
controlled  by  energy  balance. 
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IV.  Computer  Considerations 


( 


In  general,  both  of  the  programs  used  (See  Appendix  B)  do  the 
same  computations,  their  basic  difference  being  that  program  TEMPS 
plots  out  the  original  temperature  curves  for  the  air,  tank,  ground 
and  leaf  models,  and  program  DELTA  plots  the  temperature  difference 
between  the  target  (tank)  model  and  the  ground  model.  The  programs 
are  set  up  to  vary  one  parameter  over  a range  of  four  values  for 
each  run.  The  programs  must  be  changed  slightly  when  a different 
parameter  is  selected.  These  changes  amount  to  altering  the  input 
reading  cards  so  that  they  are  set  up  for  the  proper  variable  and 
putting  the  proper  title  name  into  some  of  the  graph  commands. 

Input  Data 

Input  can  be  read-in  either  by  cards  or  by  terminal.  All  input 
is  unformatted,  since  a list  directed  read,  command  is  used.  The 
organization  for  the  input  variables  is  shown  in  Appendix  B.  The 
number  of  runs  is  controlled  by  the  variable  NCRAPHS.  There  must 
be  one  complete  set  of  data  cards  (11  cards  total)  for  each  run. 

These  cards  set  up  the  initial  conditions  concerning  the  location, 
climate  and  physical  make-up  of  the  models,  such  as  relfectivity  and 
emmissivity. 

A single  run  (NGRAPHS  = 1 ) with  program  TEMPS  will  produce  a 
set  of  four  graphs  of  basic,  30-hour  temperature  profiles.  The  same 
single  run  with  program  DELTA  will  give  a single  graph  with  four 
curves  corresponding  to  the  tank -ground  ATs  of  the  graphs  in  program 
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TEMPS.  These  curves  will  start  at  the  6th  hour  because  this  minimizes 
the  initial  transient  conditions  inherent  in  the  model.  These  trans- 
ients arise  since  the  algorithm  sets  all  relevant  temperatures  at  the 
value  of  the  air  temperature  in  order  to  start  the  calculations. 


Standard  Conditions 

This  study  focuses  on  three  hypothetical  geographic  areas  with  dif- 
ferent climatic  conditions.  The  first  is  a northern  continental 
region,  the  second,  a mid-latitude  coastal  area  and  the  third  is  a mid- 
latitude  desert  location.  The  effects  of  changing  certain  parameters 
at  each  place  are  investigated  under  a variety  of  climates  designed 
to  simulate  winter  (with  and  without  snow),  spring,  and  summer.  Table  II 
summarizes  the  assumed  standard  conditions  for  each  location  and  season. 


Table  II 

Standard  Climatic  Conditions  for  Three  Model  Locations 


Location 

Latitude  Season  Temp 

Precip. 

H?0  Vapor 

(mm) 

Rel 

Hum 

{%) 

Sun 

(*) 

Dust 

Continental 

50 

/iinter(Snow) 

-1.0 

2.0 

80 

50 

0.50 

Winter (No  Snow) 

-1.0 

2.0 

80 

50 

0.50 

Spring 

7.5 

3.0 

72 

50 

1.00 

Summer 

17.0 

4.0 

66 

50 

1.00 

Coastal 

32 

iVinter(Snow) 

9.0 

1.0 

70 

50 

1.00 

Winter(No  Snow) 

9.0 

1.0 

70 

50 

1.00 

Spring 

20.0 

3.0 

67 

50 

2.00 

Summer 

27.0 

4.5 

64 

50 

_ 2.50 

Desert 

20 

rf inter 

15.0 

0.4 

40 

K 1 

Spring 

24.0 

0.6 

tq 

Summer 

22*0. 

0^2 

20 

■21 

(Ref  20i  Ref  23) 


Variables 

Not  all  of  the  31  variables  which  are  read  as  input  were  able 
to  be  examined.  Other  than  the  changes  in  location  and  climate  al- 
ready mentioned,  11  other  separate  parameters  where  run  through  the 
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( program.  Each  of  these  parameters  was  given  a number  to  simplify 

organizing  the  graph  sets  that  were  generated.  Table  III  lists  the 
parameters,  their  computer  variable  name,  their  corresponding  ident- 
ification number,  and  the  total  range  of  values  used  in  the  study. 


Table  III 

Parameter  Identification  Numbers 


I.D. 

# 

Range  of 

Values 

Parameter 

Variable  Name 

1 

-?  - +8  (C°) 

Kean  Temperature 

MEANTEMP 

2 

.16  - 7.2  (mm) 

Absolute  Humidity 

WATER 

3 

0 - 401  (cm/sec) 

Wind 

SPEED 

4 

25  - 120  (%) 

Percent  Sunshine 

CLOUD 

5 

0.08  - 1.00 

Ground  Reflectivity 

RS0IL 

6 

.475  “ 1.00 

Ground  Emlssivity 

ES0IL 

? 

.09  - .54  (cm2/8) 

Ground  Diffusivity 

DS0IL 

8 

.09  - .54 

Target  Reflectivity 

RTANK 

9 

.35  - .98 

Target  Emlssivity 

STANK 

10 

1.9  - 9.5  (cm) 

Target  Thickness 

TTANK 

11 

0 - 536  (cm/sec) 

Internal  A.  C.  Velocity 

SPEEDGI 

Plot  Identification  Numbers 

On  each  plot,  or  set  of  plots,  a specific  three -number  key  has 
been  used  to  identify  the  variable  being  examined,  the  latitude  cor- 
responding to  one  of  the  three  locations  and  the  Julian  day  that  the 
run  is  simulating.  Figure  3 shows  the  organization  of  this  system. 


f 

imete 


Parameter 
being 
varied 
(Table  III) 


PLOT  SET 

1 1 32  I IS 

'\  I 


Latitude 
of  model 
location 

50-Continental 
32 -Coastal 
20-Desert 


Julian  Day 
of  run 

1 S-W inter (Snow ) 

1 -Winter (No  Snow) 

90-Spring 

180-Summer 


Fig.  13.  Explanation  of 
Plot  Set  Numbers 
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V.  Results 


Diurnal  Temperature  Curves 

A single  run  of  program  TEMPS  results  in  four  diurnal  temperature 
plots  (one  plot  is  shown  in  Fig.  14).  All  The  temperatures  are  started 
near  the  air  temperature,  but,  by  the  fifth  hour,  each  curve  takes  on 
its  own  profile.  The  repetition  of  the  curves  at  6 and  30  hours  indi- 
cates that  the  initial  transient  condition  damps  out  within  a day.  The 
peak  tank  and  ground  temperatures  in  Fig.  14  (simulating  summer)  occur 
about  two  hours  after  the  tine  of  maximum  insolation  (at  14  hours,  versus 


Fig.  14.  Sample  Diurnal  Temperature  Curve 
for  Summer  Conditions  at  30°  Latitude. 


12  hours).  The  leaf  temperature  also  increases  above  the  air  temper- 
ature, slightly,  between  10  and  14  hours.  This  indicates  that  the 
leaf  is  transpiring  at  its  maximum  rate  and  that  its  temperature  is 
being  determined  by  energy  balance.  This  profile  parallels  nicely 
the  theoretical  curves  proposed  by  Gates  (Ref  24:82).  At  night  the 
leaf  temperature  stays  below  the  air  temperature,  an  indication  of  the 
high  amount  of  outward  radiation  ( R0  ) occuring  at  both  surfaces. 

Were  it  not  for  the  wind  warming  the  leaf  by  convection,  this  temp- 
erature drop  would  be  more  severe  (as  in  Plot  Set  3/20/280).  The 
ground;  however,  stays  slightly  warmer  due  to  the  reverse  heat  flux 
from  inside  the  earth,  and  the  tank  is  warmer  still,  mainly  due  to 
the  imposed  internal  temperature  being  constant  at  28 °G.  (The  sig- 
nificance of  this  effect  is  shown  in  Plot  Set  ll/50/l).  The  remaining 
diurnal  plots  are  located  in  Appendix  G. 

A T Curves 

Figure  15  shows  a sample  AT  plot,  where  the  target  emissivity 
has  been  varied.  From  this  graph  it  can  be  seen  that  in  general,  de- 
creasing the  target  emissivity  in  the  summer  at  northern  latitudes 
tends  to  increase  the  AT  between  the  target  and  the  ground.  Another 
characteristic  that  can  be  deduced  is  that  all  the  curves  appear  to 
have  minima  near  the  time  of  peak  heat -loading  (14.00  hours)  and  max- 
ima when  heat-loading  is  the  lowest  in  the  early  morning.  This,  too, 
is  probably  caused  by  the  preset  internal  temperature  of  the  tank  at 
28°G.  This  point  can  be  emphasized  by  noticing  the  lack  of  regular- 
ity in  the  AT  plots  for  Plot  Set  11/50/180.  Note  too  that  the  equal 
spacing  of  the  curves  in  Fig.  15  seems  to  show  a linear  relationship 
between  the  emmissivity  and  AT  for  this  computer  model.  Additional 
A T curves  are  in  Appendix  D. 
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Fig.  15.  Sample  Plot. 
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Par..,,  "ter  Sensitivity  Curves  and  Analysis 

All  of  the  parameters  that  were  varied,  except  the  mean  temp- 
erature, were  altered  around  a "standard  value"  so  that  a normalized, 
multiplicity  constant  could  be  used  to  compare  the  relative  changes 
brought  about  by  each  parameter.  In  order  to  compare  these  changes, 
the  same  two  selected  times  (14  and  24  hours)  were  used  as  data  points 
for  each  parameter.  By  graphing  the  multiplicity  value  against  the  ob- 
served AT  , the  effects  of  all  parameters  can  be  displayed  on  a single 
graph  for  each  set  of  seasons  and  latitudes.  These  plots  are  additionally 
separated  into  day  and  night  categories,  plus  the  target  and  ground/en- 
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vironraental  parameters  are  plotted  separately  to  improve  clarity. 

The  first  two  groups  of  curves  (Fig.  16  and  Fig.  17)  indicate  the 
results  that  occured  when  the  mean  temperatures  were  raised  and  lowered 
above  the  standard  values.  Figure  16  represents  the  sample  values  taken 
in  the  day.  All  of  the  curves  appear  to  have  the  same  negative  cor- 
relation, Moreover,  the  even  spacing  of  most  of  the  curves  Indicates 
that  the  same,  apparently  linear, relationship  exists  for  all  the  seasons, 
and  the  effect  of  the  season  is  just  an  upward,  or  downward  shift.  It 
is  interesting  to  note  the  fairly  large  shift  between  the  snow  and  no- 
snow lines  on  the  top  two  graphs.  It  seems  that  the  changing  of  the  re- 
flectivity of  the  ground  also  has  a strong  effect  in  shifting  the  curves. 
The  smaller  spacing  for  the  50°  latitude  example  is  apparently  due  to 
the  lower  insolation  at  that  latitude. 

The  curves  in  Fig.  17  have  all  shifted  upward  during  the  nighttime 
when  compared  to  Fig.  16.  This  is  because  the  heated  tank  exhibits  a 
stronger  signature  when  the  night  air  is  colder.  The  other  relation- 
ships mentioned  for  Fig.  16  still  seem  to  hold,  except  that  the  effects 
due  to  the  snow  are  greatly  reduced  at  night  when  there  is  no  insolation. 

Figures  18  (a)  through  18  (k)  represent  the  analysis  of  parameters 
2 through  11  (as  listed  in  Table  III) during  the  daytime.  The  graphs  are 
organized  with  a single  season  and  latitude  combination  on  each  page. 

The  left-hand  curves  are  those  that  result  from  the  ground/environment 
parameters,  and  the  right-hand  curves  are  for  those  parameters  pertaining 
to  the  target.  Comments  on  each  curve  will  be  made  on  each  individual 
page.  The  same  type  of  organization  will  apply  for  Figures  19  (a) 
through  19  (k)  which  include  the  analysis  made  on  the  parameters  during 
the  nighttime. 
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Fig.  17.  Changes  in  nighttime  AT 
caused  by  variations  in  mean  temperature 


Fig.  17(»).  Key  to  curves  in  Figs.  18(a)  to 
1900 * including  "standard"  paraaeter  values. 
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The  most  sensitive  external  parameter  is  the  wind  (b).  This  is 
primarily  due  to  the  large  temperature  difference  caused  by  the 
heated  tank  sitting  over  the  cold  ground  in  winter.  Even  a slight 
breeze  starts  to  cool  the  tank  off,  bringing  its  temperature  closer 
to  the  snow-covered  earth.  However,  as  the  wind  velocity  increases 
above  standard  its  effect  appears  to  become  less  sensitive  as  seen 
by  the  gradual  decrease  in  slope.  The  dlffusivity  has  an  opposite 
correlation,  but  a word  of  caution,  the  values  used  here  may  not  be 
appropriate  for  the  actual  dlffusivity  values  of  the  snow.  Ground 
reflectivity  has  some  variance  and  this  may  span  the  range  of  old  to 
new  snow.  Interior  convection  is  very  critical,  a stagnant  interior 
decreases  the  temperature  difference  a great  deal  even  though  the 
inside  air  may  be  warm. 


Fig.  18(b).  Parameter  Analysis}  Daytime 

Seasoni  Winter;  Latitude!  50°. 

The  relationships  in  this  figure  are  very  close  to  those  on  the 
previous  figure  except  that  the  standard  temperature  is  shifted  down 
by  about  3 degrees  when  there  is  no  snow  on  the  ground.  (This  shift 
is  similar  to  the  one  shown  in  Figures  16  and  17).  One  apparent  change 
with  the  external  parameters  is  that  the  dependence  on  the  amount  of 


sunshine  has  changed  from  slightly  positive  to  slightly  negative. 
This  is  probably  due  to  the  increased  amount  of  heat  that  the  soil 


is  picking  up,  causing  it  to  stay  warmer  for  a longer  period  of  time. 
The  ground  reflectivity  seems  to  be  even  less  sensitive  in  this 
regime.  Interior  and  exterior  convective  processes  still  dominate 
as  the  most  sensitive  parameters. 
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Fig.  18(c).  Parameter  Analysis;  Daytime 

Season i Spring;  Latitude j 50° 


The  strength  of  the  wind  parameter  has  begun  to  fall  off  as  the 
temperature  of  the  tank  gets  closer  to  the  ground  temperature.  In 
this  spring  example,  the  relative  change  brought  about  by  the  absolute 
humidity,  the  wind  and  the  amount  of  sunshine  are  practically  the  same. 
The  effect  of  the  diffusivity  is  about  the  same,  as  is  the  ground 
emimsiYlty  | however  the  ground  reflectivity  seems  to  be  getting 
more  important,  and  this  is  probably  due  to  the  increased  amount  of 
insolation  at  this  time  of  year.  The  effects  of  target  reflectivity, 
emisslvity  and  thickness  all  appear  to  have  about  the  same  relative 
merit  with  a negative  correspondence.  The  value  of  the  internal 
he&t  transfer  by  convection  still  remains  the  strongest  parameter 
of  them  all. 


Fig.  18(d),  Parameter  Analysis;  Daytime 

Season i Summer;  Latitude;  50° 

During  the  summer,  the  most  notable  change  is  that  the  importance 
of  the  internal  heat  exchange  (j)  falls  off  drastically,  this  is  primarily 
because  the  daytime  temperatures  are  becoming  closer  to  the  assumed 
28°C  of  the  inside  of  the  tank.  At  the  same  time  the  target  reflect- 
ivity is  becoming  increasing  important  with  the  higher  values  of  ' - 
Insolation.  The  wind  (b)  now  has  a positive  correlation  because  the 
temperature  of  the  tank  is  getting  warmer  than  the  air  temperature 
and  increased  values  of  the  wind  velocity  tend  to  cool  it  down  and 
drive  the  temperature  difference  toward  zero.  The  relative  strength 
of  the  other  climate  and  ground  factors  appear  to  remain  about  the 
same. 
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Fig.  18(e).  Parameter  Analysis!  Daytime 

Season «Winter(3 )Latitude i 32° 


In  the  mid-latitude  winter  with  a snow-covered  ground  the  out- 
side wind  still  is  the  most  sensitive  factor,  but  noi  nearly  as 
strong  as  was  the  case  at  50°  latitude.  The  ground  reflectivity 
has  a much  stronger  influence  at  32°  however  and  this  is  quite  likely 
due  to  the  greater  amount  of  insolation.  From  this  fact  it  is  apparent 
that  the  different  reflectivity  of  the  snow  (old  versus  new)  in  the 
lower  latitudes  may  be  able  to  change  the  values  of  Delta-T  by  as 
■uch  as  35  to  50  percent.  A similar  situation  could  occur  at  higher 
latitudes  in  early  spring.  The  other  external  factors  have  a relatively 
"standard"  appearance  as  far  as  significance  goes.  The  rate  of  convect- 
ive transfer  inside  the  tank  is  the  key  internal  factor,  although 
slightly  less  important  than  at  50°  with  a snow  covered  ground.  This 
is  due  to  the  warmer  air  temperature  most  likely. 
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rig.  Ib^f ).  Parameter  Analysis;  Laytime 

Season;  Winter;  Latitude;  32°. 

Without  snow,  the  ground  absorbs  more  solar  energy  and  hence 
is  warmer  during  the  daytime,  thus  the  Lelta-T  is  lower  in  this  in- 
stance and  the  effects  of  the  wind  (b)  and  the  ground  reflectivity  (d) 
have  becone  less  significant.  3ecause  the  ground  is  absorbing 

more  energy* the  diffusivity  and  emmissivity  have  become  more  sensitive, 
since  they  are  apparently  acting  as  regulators  for  the  energy  flow 
into  and  out  of  the  ground.  Very  little  change  has  occured  with  the 
internal  parameters,  in  fact. the  convective  heat  transfer  appears  to 
be  as  important  as  it  was  when  there  was  snow  on  the  ground. 


Fig.  18(g).  Parameter  Analysis;  Jayuime 

Season:  Spring;  Latitude:  32° 


As  th®  year  goes  on  and  the  air  temperature  warms  up,  the  cooling 
effect  of  the  wind  on  the  tank  becomes  a warming  effect,  as  shown  in  the 
change  in  slope  of  curve  (b)  compared  to  the  winter.  The  negative 
slope  for  the  amount  of  sunshine  (c)  indicates  that  with  increased 
insolation,  the  ground  heats  up  more  than  the  tank.  A similar  effect 
occurs  when  the  absolute  humidity  (a)  is  increased.  The  increased 
spring-time  air  temperature  causes  the  internal  heat  convection  tern 
to  become  much  less  effective,  since  the  interior  tank  temperature 
is  closer  to  the  air  temperature.  The  emissivity  of  the  tank  is 
less  important  now,  due  to  a higher  percentage  of  insolation  versus 
atmospheric  radiation  (as  compared  to  the  winter  heat-load). 


With  maximum  insolation,  the  factors  that  tend  to  cause  heat- 
gain  in  the  ground  are  showing  a strong  positive  correlation,  this 
includes  the  enisslvity.  , reflectivity  and  diffusivity.  Additionally, 
as  the  wind  velocity  decreases  outside, the  ground  heats  up  more  than 
the  tank  causing  the  value  of  the  Delta -T  to  become  more  negative. 

The  negative  slope  on  the  internal  A.G.  velocity  term  shows  that 
during  the  middle  of  the  day  the  inside  of  the  tank  is  cooler  than  the 
ground  and  cooling  the  inside  shows  up  as  causing  the  Delta-T  value  to 
become  more  negative.  An  explanation  of  why  increasing  the  tank 
thickness  causes  the  Delta-T  to  become  more  negative  (the  tank  skin 
la  cooler)  may  be  that  due  to  the  increased  heat  storage  capability, 

the  peak  tank  temperature  shifts  to  a time  that  increases  the  temp- 
erature difference  at  the  14  hour  point. 
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The  first  noticeable  part  of  the  winter  curves  in  the  desert  is 
that  the  wind  has  very  little  effect  upon  the  thermal  signature.  Cn 
the  contrary,  the  soil  reflectivity  is  of  much  greater  importance, 
now,  compared  to  the  other  two  latitudes  in  winter.  This  is  probably 
due  to  the  greater  amount  of  insolation  for  the  time  of  year.  More 
evidence  of  the  relatively  strong  Insolation  is  the  Importance  of  the 
ground's  diffusivity  and  emmissivity.  The  negative  slope  on  the  sun 
strength  curve  (c)  shows  that  as  the  insolation  is  increased  the 
heating  of  the  soil  is  greater  than  the  tank,  causing  a decrease  in 
the  Delta-T  despite  the  fact  that  the  reflectivity  of  the  soil  is 
higher  in  this  model  than  the  others  (assuming  that  it  is  desert 

sand).  With  higher  air  temperature  and  insolation,  the  tank  emissivity, 
thickness,  and  reflectivity  ere  more  important,  while  the  internal  con- 
vection is  less  important  than  at  the  other  latitudes  studied. 
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As  the  intensity  of  the  insolation  increases  from  winter  to 
spring,  the  importance  of  the  soil  reflectivity  increases  also.  A 
strong  positive  slope  on  the  wind  term  indicates  that  cooling  of  the 
sand  by  a breeze  can  change  the  value  of  the  temperature  difference 
by  a great  deal.  The  other  external  parameters  are  similar  to  the 
winter  case.  The  action  of  the  internal  convection  term  now  tends 
to  cool  the  outside  of  the  tank  and  thus  when  it  is  increased  the 
magnitude  of  Delta-T  also  decreases.  The  slope  of  the  target  re- 
flectivity has  become  slightly  more  negative,  which  is  Indicative  of 
the  increase  in  solar  insolation  over  the  winter. 
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Under  maximum  insolation  the  extreme  sensitivity  of  the  ground 
reflectivity  is  readily  apparent*  a 100#  change  in  the  reflectivity  cam 
mean  a 200#  change  in  the  Delta-T  value.  At  the  same  time,  the  cooling 
effect  of  the  wind  on  the  soil  has  also  become  a bit  more  important, 
while  the  relative  slopes  of  the  other  4 external  factors  have  remained 
about  the  same.  The  cooling  effect  of  the  internal  air  conditioning 
is  at  a maximum  and  the  difference  between  having  it  off  and  having 
it  on  can  mean  a 170#  difference  in  the  Delta-T  value.  The 
importance  of  the  reflectivity  of  the  tank  has  also  become  more 
important  with  the  increased  heat-load  of  the  summer. 

• j ; 
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Fig.  19(a).  Parameter  Analysis;  Nighttime 

Season:  'rfinter(Snow) ; Latitude;  50°. 


The  two  most  important  external  factors  at  night  appear  to  be 
the  wind  and  the  emissivity  of  the  snow,  the  former  is  because  of 
the  cooling  of  the  tank  that  will  occur  as  the  wind  velocity  increases, 
the  latter  is  due  to  the  fact  that  thermal  radiation  from  the  sky  and 
snow  comprise  the  major  energy  balance  at  the  snow's  surface.  Except 
for  the  target  reflectivity,  the  internal  terms  appear  to  have  more 
significance  at  night.  Again,  because  of  the  importance  of  thermal 
radiation  at  night,  the  eaissivity  plays  an  important  role  in  the 
target  signature.  The  positive  correlation  with  thickness  is  pre- 
sumably an  indication  of  the  gradual  release  of  the  energy  stored 
up  within  the  skin  of  the  tank. 
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Fig.  1 9(t> ) . Parameter  Analysis;  Nighttime 
Season;  Winter;  Latitude;  50. 

There  is  no  apparent  change  in  the  energy  balance  relationships 
using  a scenario  without  snow  under  the  same  conditions  as  in  Fig.  19(a). 
The  graphs  do  look  a bit  different,  but  this  is  due  primarily  to  a change 
of  scale  between  the  two. 


Fig.  19(c).  Parameter  Analysis;  Nighttime 
Season*  Spring;  Latitude*  50° 


With  Spring,  the  importance  of  the  wind  begins  to  decrease  as  the 
air  temperature  and  tank  temperature  begin  to  coincide.  The  ground 
•aissivity  has  about  the  same  effect  as  in  the  winter,  but  the  dif- 
fusivity  appears  to  be  even  less  sensitive, even  though  the  day  time 
heat -loading  is  more  at  this  time.  It  is  interesting  to  note  that 
the  intensity  of  the  sun  seems  to  have  some  effect,  only  if  it  is  re- 
duced below  the  normal  value.  It  seems  that  with  less  sun  strength, 
the  ground  collects  less  energy,  is  cooler  at  night, and  the  Delta-T 
is  thus  increased.  The  only  significant  change  with  the  target  par- 
ameters is  that  the  strong  changes  due  to  the  convective  transfer 
are  beginning  to  ebb,  primarily  due  to  the  smaller  difference  between 
the  air  temperature  and  the  target. 
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Fig.  19(d).  Parameter  Analysis;  Nighttime 
Season : Summer 5 Latitude;  50°. 

The  effects  of  the  ground's  enissivity  and  reflectivity  do  not  seem 
tc  have  changed  much, and  the  convection  term  has  nearly  zeroed  out  in 
the  summer  night.  There  is  still  a slight  negative  correlation  between 
the  variation  of  the  sun's  intensity,  but  not  as  much  as  in  the  spring. 
The  biggest  change  is  due  to  the  diffusivity  and  the  absolute  humidity. 
Apparently,  with  a lower  diffusivity  there  is  less  energy  stored  deep 
in  the  soil  to  allow  it  to  be  warmer  at  night  (the  heat  radiates  away 
faster  causing  the  nighttime  minimum  to  occur  at  a different  time).  The 
decrease  in  the  absolute  humidity  will  tend  to  decrease  the  long-wave 
radiation  at  night*  but  during  the  day  it  would  have  caused  an  increase 
in  the  insolation,  so  that  the  net  effect  is  that  the  tank  receives  more 
energy  in  the  day,  than  the  soil  receives  at  night,  thus  the  Delta -T  is 
increased.  The  convection  inside  the  tank  causes  little  change  now, 
but  the  thickness  becomes  more  Important,  as  there  is  more  solar  energy 
available  to  store  up  during  the  day. 
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Fig.  19(e).  Parameter  Analysis}  Nighttime 

Season j Winter  (Snow);  Latitude :32  . 


The  ground  eaissivity  and  wind  curves  are  similar  to  those  for 
the  spring  at  latitude  50°,  and  the  absolute  humidity  curve  has  a 
small  negative  slope  as  was  the  earlier  case.  There  are  only  very 
slight  effects  due  to  changes  in  the  diffusivlty  and  the  amount  of 
sun  insolation.  There  is  an  interesting  dip  in  the  reflectivity 
curve  for  values  less  than  about  .6  of  the  normal  value j however, 
this  is  even  a bit  low  for  a reflectivity  for  old  snow.  The  graph 
of  the  internal  parameters  is  similar  to  the  other  night  plots  in 
the  winter,  with  strong  influences  due  to  convection,  thickness  and 
eaissivity  . 
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Fig.  19(f).  Parameter  Analysis}  Nighttime 
Seasoni  Winter;  Latitude i 32°. 

The  external  factors  of  diffusivity,  sunshine,  and  reflectivity 
have  become  even  more  neutral  compared  to  the  winter  simulation  with 
snow.  The  effect  of  the  absolute  humidity  appears  to  be  stronger 
at  this  time,  more  like  it  was  for  the  summer  in  the  50°  graphs,  this 
may  be  due  to  the  increased  amount  of  insolation  at  this  latitude. 
There  has  been  very  little , if  any,  change  in  the  internal  parameter 
curves,  with  the  internal  convection  being  the  most  sensitive  still 
(due  mostly  to  the  cold  night  temperatures). 
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Pig.  19^g).  Parameter  Analysis;  Nighttime 
Season:  Spring;  Latitude:  32°. 

The  major  change  in  these  curves  is  the  loss  of  importance  of 
the  convection  terns.  This  is  most  likely  due  to  the  warmer  air 
temperatures,  and  a subsequent  drop  in  the  temperature  difference 
between  the  tank  and  the  air  inside  and  outside  the  surface.  The 
only  other  noticeable  change  is  the  increase  in  sensitivity  of  the 
dlffusivity  of  the  soil.  This  is  probably  caused  by  the  higher 
Insolation,  and  thus  a higher  heat-load  at  the  surface  of  the  ground. 
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Fig.  19(h).  Parameter  Analysis:  nighttime 
Season:  Summer}  Latitude:  32°. 

With  the  warm  summer  evenings  the  slopes  of  the  convection  terms 
continue  to  change  as  the  night  air  has  a tendency  to  warm  the  tank 
and  the  air  conditioning  within  cools  it.  The  ground  diffusivity 
parameter  still  becomes  more  effective  as  the  daily  insolation  now 
reaches  its  peak  so  that  the  sun  is  up  longer  and  more  energy  gets 
stored  in  the  soil  to  radiate  (and  convect)  away  at  night.  The  slopes 
of  the  smissivities  of  both  the  tank  and  the  ground  appear  to  be  a 
little  less  steep.  This  is  probably  due  to  the  warmer  air  temperatures 
causing  more  atmospheric  thermal  radiation  so  that  the  net  cooling  effect 
due  to  radiation  is  less. 

< 
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Fig.  19(l).  Parameter  Analysis:  nighttime 
Season:  Winter;  Latitude:  20°. 

The  angles  of  the  desert  curves  are  fairly  similar  to  those  at 
32°  in  the  winter,  except  that  they  are  shifted  down  somewhat  (the 
air  temperature  is  warmer)  and  that  the  convective  heat  terms  are  not 
nearly  so  dominant,  again  because  of  the  warmer  air  temperature . There 
is  a slightly  more  downward  slope  to  both  the  sunshine  (c)  and  the 
dlffusivity  (f)  curves  than  there  was  in  the  previous  latitude  in  the 
winter.  This  may  be  due  to  a higher  amount  of  insolation  during  the 
day  accentuating  the  effects  of  these  terms  during  the  nighttime.  This 
may  even  affect  the  reflectivity  curve  of  the  ground  as  well. 

The  target  reflectivity  still  indicates  an  almost  flat  response. 
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Fig.  1 9( j ) . Parameter  Analysis:  nighttime 

Season:  Spring}  Latitude:  20°. 

i 

Once  again  the  convection  terms  begin  to  swing  to  different 
slopes  as  the  air  temperature  warms  up,  and  again,  as  before,  the  soil 
diffusivity  slope  increases  negatively.  There  also  appears  to 
be  a very  slight  change  in  the  sensitivity  of  the  ground  reflectivity 
but  not  very  much.  Except  for  the  convection  term,  the  internal 

factors  are  almost  totally  unaltered  in  the  change  of  season  from 
winter  to  spring.  The  upward  displacement  of  the  absolute  humidity 
terms  is  possibly  due  to  a slight  change  in  the  computer  program  when 
that  date,  set  was  run.  There  was  insufficient  time  to  go  back  and 
verify  the  data.  At  any  rate  the  size  of  a vertical  shift  in  the 
curve  is  rather  arbitrary  since  the  changes  in  slope  are  what  is 
most  important  for  this  analysis. 

1 
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Fig.  19(kJ.  Parameter  Analysis : nighttime 
Season:  Summer;  Latitude:  20°. 

Because  of  the  large  difference  between  the  tank  temperature  and 
the  air  during  a desert  summer  night  the  convective  terms  dominate  as 
the  most  important  factors.  Were  the  basic  temperature  inside  the  tank 
changed,  a subsequent  change  in  the  convective  terms  will  show  up  8s 
well.  As  was  also  seen  in  the  set  of  curves  for  the  summer  at  32° 
latitude,  the  diffusivity  of  the  soil,  and  the  eniBsivities  of  the 
soil  and  tank  have  flattened  out  somewhat,  possibly  another  result  of 
the  increased  air  temperature.  No  other  significant  changes  seem 

to  have  occured. 
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VI.  Conclusions  and  Recommendations 


Conclusions 

Table  IV  Is  a summary  of  ranking  the  parameters  using  the  curves  in 
Figs  18(a)  to  19(k).  The  letters  refer  to  the  parameters  (see  Fig. 
!?(*)) t and  ar®  listed  from  high  to  low  Importance  (left  to  right) 


Table  IV 


Summary  of  Parameter  Strengths  for  Eleven 
Hypothetical,  Geographic  and  Climatic  Conditions 


Condition 

Latitude 

Day 

Night 

Winter 

(Snow) 

50°N 

32°N 

j *b,f ,e,d,h,i,a,g,c 
j*b,d,f ,h,i,e,c,g,a 

j»b,e,h,i,f ,a,d,c,g 
b» j,e,i,h,a,d,f ,c,g 

Winter 
(No  Snow) 

50°N 

32°N 

20°N 

j»b,f ,e,h,i,d,a,c,g 
J»b,e,f ,h,i,g,d,c,a 
di®,f,g,h,i,j,c,a,b 

j»b,e,h,i,f ,a,d,c,g 
j,b,e,h,l,a,f ,d,c,g 
e,i, j,h,b,c,f ,a,d,g 

Spring 

50°N 

32°N 

20°N 

j,f,e,b,g,h,c,i,a,d 
f,e,d,g,a,c,h,i,b,j 
d|b,f ,e,g, j,h,l,e,a 

j»b,h,e,l,a,c,d,f ,g 
®#i*a»h,f ,b,c,d, j,g 
e,i,f,h,j,b,d,a,c,g 

Summer 

50°N 

32°N 

20°N 

f »®*g*d»h,i,a,c,b, j 
f,d,b,e,j,g,i,h,a,c 
d.j.g.b.e.f .h.c.l.a 

®,i,a,h,f ,c,d, j,b,g 
j,e,a,f,i,b,h,g,d,c 

Jj®.rli^|TAhlaJd,cJg 

If  values  of  1 (most  Important)  to  10  (least  important)  are  as- 
signed to  the  ranking  of  each  scenario,  and  if  these  values  are  sum- 
med for  each  parameter,  an  overall  parameter  strength  can  be  deter- 
mined (the  lowest  sun  representing  the  highest  overall  importance). 

The  day  and  night  totals  for  the  parameters  are  shown  In  Table  V. 

Although  ground  emisslvity  is  the  most  important  overall  factor, 
its  value  changes  very  little  for  most  natural  objects  in  the  8 to  12 
micron  region  (See  Appendix  E).  On  the  other  hand,  the  diffusivity 
varies  greatly  for  different  soils,  and  the  wind  velocity  changes  con- 
tinually, so  these  parameters  can  be  the  sources  of  significant  changes 
in  Delta-T.  The  low  insolation  rating  is  surprising,  but  may  actually 
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Table  V 


Overall  Day  and  Night 

Parameter  Strength  Rating 

Day  Night 

(i.D. ) Parameter  Name  (Sum)  (i.D. ) Parameter  Name  (Sum) 


(fj 

Ground  Diffusivity 

31 

Ground  Emissivity 

24 

ie( 

Ground  Emissivity 

40 

j) 

Internal  A.C.  Velocity  33 

(d 

Ground  Reflectivity 

4 5 

1 

Target  Thickness 

40 

(j, 

Internal  A.C.  Velocity  45 

,b) 

Wind 

45 

Wind 

49 

h) 

Target  Emissivity 

49 

(g 

Target  Reflectivity 

65 

[a) 

Absolute  Humidity 

64 

(h) 

Target  Emissivity 

66 

,f) 

Ground  Diffusivity 

65 

u; 

Target  Thickness 

77 

(d) 

Ground  Reflectivity 

87 

(0) 

Solar  Insolation 

92 

lc) 

Solar  Insolation 

90 

(a, 

Absolute  Humidity 

95 

ig) 

Target  Reflectivity 

108 

be  caused  by  the  small  difference  (.1)  between  the  reflectivities 
of  the  standard  target  and  background.  The  insonsitivity  of  the 
absolute  humidity  may  be  due  to  the  assumption  of  drier  than  normal 
standard  climates.  The  data  in  Table  V should  be  used  cautiously, 
because  many  seasonal  climate  characteristics  have  been  encompassed 
into  one  figure  of  merit. 


Rec  ommendatlons 

To  improve  the  reliability  of  the  program,  the  basic  weather  models 
should  be  made  more  analytical  to  allow  better  adjustments  for  meteor- 
ological conditions.  The  test  models  (the  ground,  tank,  and  trees) 
should  be  checked  to  ensure  that  the  heat  transfer  mechanisms  are 
represented  in  the  proper  magnitude. 

An  attempt  should  be  made  to  validate  the  computer  model  by  running 
It  with  actual  real-time  data,  suoh  as  temperature,  wind,  cloud-cover, 
and  insolation.  A tine  increment  of  15  minutes  is  recommended  for  such 
a trial,  since  it  seems  to  work  well  for  the  heat  transfer  models  in 
this  study.  Longer  increments  would  cause  poorer  resolution,  and  short- 
er increments  lead  to  problems  where  time  independence  is  assumed. 
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Appendix  A 


Comparison  of  Ground  Conduction  Model  to 
Analytical  Solution  Using  Sinusoidal  Input 

Fig.  20  shows  the  analytical  solution  to  ground  temperature  for 
a semi -infinite  slab  as  computed  from  Carslaw  and  Jaeger  (Ref  21:76). 
Fig.  21  shows  the  solution  using  Eq.  (14)  and  a time  increment  of  .25 
and  Fig.  22  shows  the  Eq.  (14)  solution  with  a time  increment  of  .02. 
Additionally  Fig.  23  gives  an  example  of  the  TEMPS  program  plot  with 
a .25  hour  time  increment  and  Fig.  24  gives  an  example  of  another 
TEMPS  plot  with  a .15  hour  time  increment.  The  greatest  difference 
between  Fig.  23  and  Fig.  24  is  the  maximum  of  the  ground  temperature  j 
it  is  equal  to  30.89  °G  when  DT-.25  and  it  equals  31.07  when  DT“.15. 
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Fig.  21  Eq.  (14)  solution  to  semi -infinite  slab 
problem,  with  DT-.25. 
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Fig.  22  Eq.  (14)  solution  to  semi-infinite  slab 
problem,  with  DT-.02. 
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Appendix  B 

Computer  Programs  TEMPS  and  DELTA 


Input  Data  t 
Card  Symbol 

1 NGRAPHS 

DMIN 

LL 


LAT 

DAY 

SPEED 

CTEMP1 

NCLOUD 


CLOUD 


WATER 


RH 


Description 

Determines  how  many  sets  of  plots  are  to  be  plotted. 
Sets  the  minimum  value  for  the  graphs  in  program 
DELTA  only.  Is  a dummy  variable  in  program  TEMPS. 
An  integer  which  determines  whether  the  letter  "S" 
is  printed  after  some  of  the  plotting  information 
to  indicate  that  the  data  Is  set  for  snow  on  the 
ground.  (Data  is  read  in  by  the  operator  though) 

The  latitude  in  degrees  of  the  model  location. 

The  day  of  the  run  (Julian  calendar). 

Velocity  of  the  wind  in  cm/sec 
DC  term  in  air  temperature  calculation  °C 
An  integer  that  determines  the  type  of  clouds 
present  it  corresponds  to  the  grouping  listed 
in  Table  I and  sets  the  appropriate  values  for 
the  particular  type  of  cloud. 

A real  number  from  0.0  to  1.0  to  indicate  what 
fraction  of  the  insolation  actually  reaches  the 
earth,  a type  of  filter  factor. 

Real  number  which  indicates  the  total  precip- 
1 table  water  vapor  in  a column  one  meter  square 
reaching  through  the  atmosphere. 

Relative  Humidity  in  a number  from  0 to  100. 
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Card 

Symbol 

Description 

2 

PRESS 

Atmoshperic  pressure  in  millibars. 

DUST 

A measure  of  the  particulate  matter  in  the  air  in 

particles/cc. 

3 

DESOIL 

The  density  of  the  soil  in  gm/cc. 

CPSOIL 

Specific  heat  of  the  soil  in  Cal/gm-°C. 

DSQIL 

Diffusivity  of  soil  in  cm^/min. 

ESOIL 

Emmissivity  of  the  soil. 

RSOIL 

Reflectivity  of  the  soil. 

KSOIL 

Thermal  conductivity  of  the  soil  in  Cal/min-cm-°C. 

WET 

The  amount  of  moisture  in  the  soil  from  0.0  to  1.0 

k 

DETANK 

The  density  of  the  skin  of  the  tank  in  gm/cc. 

CPTANK 

Specific  heat  of  the  tank  material  in  Cal/gm-°C. 

ETANK 

Emmissivity  of  the  tank. 

RTANK 

Reflectivity  of  the  tank. 

SPEEDGI 

Velocity  of  air  movement  inside  tank  in  cm/sec. 

TEMFGI 

Internal  temperature  of  the  tank  °C 

TTANK 

Thickness  of  the  tank  in  cm. 

5 

DELEAP 

The  density  of  the  leaf  in  gm/cc. 

CPLEAF 

Specific  heat  of  the  leaf  in  Cal/gm-°C. 

DIMEN 

Dimension  of  leaf  along  the  flow  of  air  in  cm. 

TRANS 

The  maximum  rate  of  transpiration  of  the  leaf  in 

grams  per  min. 

ELEAF 

Emmissivity  of  the  leaf. 

A LEAF 

Absorptivity  of  the  leaf. 

TLEAF 

Thickness  of  the  leaf  in  cm. 
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Symbol  Description 

DT  Length  of  the  time  increment  in  hours. 

HRS  Length  of  desired  plot  in  hours, 

PLOTNO  Integer  used  to  identify  the  set  of  plots  with  a 
certain  parameter  that  is  varied  (See  Fig.  13) 

DSOIL  Any  parameter  that  is  being  varied  (in  this  case  it 
is  the  diffusivity  of  the  soil).  This  can  be  any  of 
the  desired  values  that  the  operator  choses. 

XXX  The  minimum  value  of  the  parameter  to  be  varied. 

This  value  is  printed  in  a block  on  the  plot  set 
that  is  produced  by  program  TEMPS. 

YYY  The  maximum  value  of  the  parameter  to  be  varied. 

DSOIL  Second  value  of  parameter  to  be  used  in  second  time 
through  the  program. 

DSOIL  Third  value  of  parameter  for  third  run. 

DSOIL  Fourth  value  of  parameter  for  fourth  run. 

The  sequence  of  cards  2 through  10  is  repeated  as  many  times 
as  the  number  NGRAPHS  has  been  set  to. 

The  remainder  of  this  appendix  contains  a listing  of  the  Fortran  pro- 
gram, TEMPS  including  all  subroutines,  then  the  Fortran  program  DELTA, 
excluding  subroutines  (they  are  identical  to  those  used  in  TEMPS)  and 
then  a sample  output  listing  from  program  DELTA. 
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*******»*****»********.»**********♦** 

INFOPMATION  FCR  PLOT  NUMBER  - 2. 
ENVIRONMENTAL  rACTORSt 


LATITUDE*  20. 

OAY=  so. 

WIND  <°FEO=  223. 
MEAN  T EMP=  24. 
CLO'JO  TYPE  = 1 
PFRCFNT  SJN=  .T 
HATFP  VAPDRs  .6 
HUHIOITYs  ?0. 
PRESSURES  10  0 0. 
OUST=  .25 


GROUND  °ARAMETERS 


DFNSITYs  2. 

SPECIFIC  HEAT=  .44 
DIFEUSIVITYs  .3 
EMISSlVITYs  .<35 

RFFLFCTIVITYs  .3 

CONOUCTTVITYs  . ?i 
MOISTURE  CONTENT-  .15 

TANK  PARAMETERS 


OENSITYs  B . 

SPECIFIC  HEAT=  .11 
EMISSlVITYs  .v 
REFLECTlVITYs  .3 

AIP  CONDITIONING  VELOCITYs  44T. 
INTERNAL  T FMPER \TURE=  25. 
THICKNFSS=  7.5? 

LEAF  PARAMETERS 


OENSITYs  .5 
SPECIFIC  HFATs  .as 
SI 7E=  2.54 

TRANSPIRATION  PATEs  .0003 
EMISSlVITYs  .95 
ABSORPTION:  .7 
THlNCKFqs  =.i 

***«*#»*#***» *» ************************ 


TOTAL  INSOLATION  WASt  473.2973654092  CAL/3Q.CM. 
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DATA  FOR  ©LOT  S"T  2./20./  90. 


PEAK  TEMPS  ANT  "'IMES-I  TERATION  #1 
********************************** 


AIR 

34.  68  0*0^78449 

16.5 

GROUND 

33.  87  0 ’05  21 637 

14. 

TANK 

34.919  1131613 

14.5 

LEAF 

34. 59 93 0426596 

16. 

********************************** 

* 1400 

DELTA -Tt 

.84  ♦ 

♦ 240  0 

DELTA-T*  5 

.10  * 

********  ********  ******************* 

PEAK  TFMPS  AND  timfS-ITERATION  *2 

********************************** 

AIR 

34.  680  10’78449 

16.5 

GROUND 

33.  766  18620909 

14. 

TANK 

34.  84 2 (4163599 

14.5 

LEAF 

34.  58e  <04  26596 

15. 

* 1400  OELTA-T!  .87  * 

* 2400  DELTA  -T!  6.10  * 

PFAK  TFMPS  ANT  "IMES-ITERATtON  #3 


AIR 

34.  69  040',78449 

16. 

5 

GROUND 

33.  684e4322091 

14. 

TANK 

34.  ■”7986380326 

14. 

5 

L EAr 

34.  58580426596 

16. 

************* >*****+************** 

* 1400  DFLTA  -Tt  .89  * 

* ?UOO  OFLTA-Ti  5.11  * 

*********************************** 

PFAK  Tfmps  AND  ~I MES -I TERA T ION  #4 


AIR  34. 68040’78449 

16.5 

GROUND  33.61412670207 

14. 

TANK  34.  721  T4D13484 

14.5 

LEAF  34.58590*26596 

16. 

♦ 1400  DELTA-T! 

.90 

* 

* 2400  OELTA-Tt  6 

.11 

* 

i 
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Appendix  C 


Plots  from  Program  TEMPS  for  all  Parameters 
Figs.  25(a)-25(k)  are  plots  of  the  parameter  Mean  Temperature. 
Figs  26(a)-26(k)  are  plots  of  the  parameter  Absolute  Humidity.  Figs. 
27(a)-2?(k)  are  plots  of  the  parameter  Wind.  Figs  28(a)-28(k)  are 

plots  of  the  parameter  Sun  Strength.  Figs.  29(a)-29(k)  are  plots 

/ 

of  the  parameter  Ground  Reflectivity.  Figs  30(a)-30(k)  are  plots  of 
the  parameter  Ground  Emmissivity.  Figs  3l(a)-3l(k)  are  plots  of  the 
parameter  Ground  Diffusivity.  Figs.  32(a)-32(k)  are  plots  of  the 
parameter  Target  Reflectivity.  Figs  33(a)-33(k)  are  plots  of  the 
parameter  Target  Emmissivity,  Figs  3^(a)-3^(k)  are  plots  of  the  par- 
ameter Target  Thickness.  Figs.  35(a)-35(k)  are  plots  of  the  parameter 
Internal  Air  Conditioning  Velocity. 


TCnPEMTUHE  fOCO  C» 


Fig.  26(e).  Basic  Temperature  Plots  of  Absolute  Humidity 
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Fig.  27(f).  Basic  Temperature  Plots  of  Wind 
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Fig.  28(d).  Basic  Temperature  Plots  of  Sun  Strength 
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TenrcMTURE  coeo  c> 


Fig.  28 (i).  Basic  Temperature  Plots  of  Sun  Strength 
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Fig.  29(a).  Basic  Temperature  Plots  of  Ground  Reflectivity 
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TCRrCMrTUM 


TcnrCMTUM 


Fig.  29(h).  Basic  Temperature  Plots  of  Ground  Reflectivity 
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UfirCMTUM  COCO  C)  TOIPCMTUNE  COCO  C) 


Fig.  29 (k).  Basic  Temperature  Plots  of  Ground  Reflectivity 
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Fig.  30(e).  Basic  Temperature  Plots  of  Ground  Emmisslvity 
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Fip.  30(f),  Basic  Temperature  Plots  of  Ground  Emmissivity 
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Fig.  30(h).  Basic  Temperature  Plots  of  Ground  Emraissivity 


Fig.  30(i).  Basic  Temperature  Plots  of  Ground  Emraissivity 


Fig,  3l(*).  Basic  Temperature  Plots  of  Ground  Diffuslvity 
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TEfirEMTUHE  (OEO  CJ 


TCftPElWTURE 
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Fig.  31 (k).  Basic  Temperature  Plots  of  Ground  Diffusivlty 
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«■  • **'.  — perat  .re  flat*  of  Target  Reflectivity 
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Fig.  32(c).  Basic  Temperature  Plots  of  Target  Reflectivity 
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Fig.  32(f).  Basic  Temperature  Plots  of  Target  Reflectivity 


170 


TcnrcMTbM 


Fig.  32(h).  Basic  Temperature  Plots  of  Target  Reflectivity 
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Tig.  32(l).  Basic  Temperature  Plots  of  Target  Reflectivity 
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Fig.  32 (k).  Basic  Temperature  Plots  of  Target  Reflectivity 
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Fig.  33(d)‘  Basic  Temperature  Plots  of  Target  iinunissivity 
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Fig*  33(1).  Basic  Temperature  Plots  of  Target  Emmlssivity 
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Fig.  3^(d).  Basic  Temperature  Plots  of  Target  Thickness 
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Fig.  y+(g).  Basic  Temperature  Plots  of  Target  Thickness 
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Fig.  3^(h).  Basic  Temperature  Plots  of  Target  Thickness 


35(c).  Basic  Temperature  Plots  of  A.C.  Velocity 
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fig*  35(d).  Basic  Temperature  Plots  of  A.C.  Velocity 
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Fig.  35(f).  Easic  Temperature  Plots  of  A.C.  Velocity 
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Fig.  35(k).  Basic  Temperature  Plots  of  A.C,  Velocity 
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Appendix  D 

Plots  from  Program  DELTA  for  all  Parameters 
K 36(a)-36(d)  are  plots  of  the  parameter  Mean  Temperature. 
Figs.  37(a)-3?(d)  are  plots  of  the  parameter  Absolute  Humidity.  Figs. 
38(a)-38(d)  are  plots  of  the  parameter  Wind.  Figs.  39(a)-39(d)  are 
plots  of  the  parameter  Sun  Strength.  Figs.  40(a)-40(d)  are  plots  of 
the  parameter  Ground  Reflectivity.  Figs  41 (a)-4l (d)  are  plots  of  the 
parameter  Ground  Emmissivity.  Figs  42(a)-42(d)  are  plots  of  the  par- 
ameter Ground  Diffusivity.  Figs.  43(a)-43(d)  are  plots  of  the  par- 
ameter Target  Reflectivity.  Figs  44(a)-44(d)  are  plots  of  the  para- 
meter Target  Emmissivity.  F'igs45(a) -45(d)  are  plots  of  the  parameter 
Target  Thickness.  Figs  46(a)-46(d)  are  plots  of  the  parameter  In- 
ternal Air  Conditioning  Velocity. 
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Fig.  36(c).  Delta-T  Plots  for  Meantemp 
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Fig.  36(d).  Delta-T  Plots  for  Meantemp 
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Fig.  37(a).  Delta -T  Plots  for  Absolute  Humidity 
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Fig.  41(c).  Delta-T  Plots  for  Ground  Emmissivity 
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Fig.  43(d).  Delta-T  Plots  for  Target  Reflectivity 
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Fig.  46(d).  Delta -T  Plots  for  A.C.  Velocity 
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Appendix  K 


Staple 


Duna  Physical  Characteristic! 


Table  VI  ia  a Hating  of  aoae  of  the  physical  characteristics 
of  the  targets  and  backgrounds  used  In  this  report*  they  are  froa 
a variety  of  references*  listed  in  the  Bibliography.  Figure  4?» 
taken  froa  he Glatc hey  (Ref  29«75)  shows  the  variation  of  the  reflec 
tivitles  of  aoae  natural  objects  in  the  range  of  .4  to  14  ai crons. 


Table  VII 


Sample  Target/Background  Charaoterl sties 
Characteristic  Substance  Value  Reference 


Reflectivity 


(Albedo) 


Soil (Vet) 

Soil (Dry) 

Saudi Vet) 

Sand (Dry) 
•Terrain" 
Oxidised  Iron 
Brown  Paint 
Olive  Drabe  Paint 
Oak  Leaf 
Haple  Leaf 
Grass 


.05 

.16 

.008 

.087 

.2 

.1 

.01 

.06  - .07 
.1  - .45 
.1  - .4 
.25 


Ref  14sl4 

m 

m 


Ref  14i8 

Ref  27 
« 


Ref  26tl4 


Diffusivity 


Soil (Dry) 

Soil (Vet) 
Soil(Avg) 

Soil (Sandy, Dry) 
Soil (8*  H20) 
Soil (Damp) 
Sandstone 

n 


Lieeetone 


Rock(Avg) 

Rock 

Snow 

Snow (Fresh) 

Vater 

Ioe 


.01  ftz/hr(.l5  cm2/min)  Ref  16*594 
.0031  cmZ/sec(.186  cm  /min)  Ref  28*143 
.0046  cm2/sec(.276  cm2/min)  Ref  21*497 
.002  cm2/sec(.12  cm2/min)  " 

.0033  cm2/sec(,198  cmr/min)  " 

.005  cm2 /sec (.3  cm2/min)  Ref  28*143 
.0113  cm2/sec(.678  cm2/min)  " 

.011  cm2/sec(.66  cm2/min)  Ref  21*497 
.045  ft2/hr(.7  cm2/min)  Ref  15*508 
.017  ft2/hr(.26  cm2 /min)  Ref  16*594 
.007  cm2/sec(.42  cm2 /min)  Ref  21  *497 
.0081  cm2/sec(.49  cm2/ min)  Ref  28*143 
.0118  cm2/sec(.71  cm-Zmin)  Ref  21*497 
.0118  cm2/ sec (.71  cm2/min)  Ref  28*143 
.005  cm2/sec ( . 3 cm2/min)  Ref  21*497 
.0033  cm2/sec(.198  cmz/min)  Ref  28*143 
.00144  cm2/sec  ( . 084  cm2/mln)  Ref  21*497 
.0115  cm2/sec(.68  cm2/min)  " 
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